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MASTER  DATA  BASE  FOR  OPTICAL  TURBULENCE 
RESEARCH  IN  SUPPORT  OF  AIRBORNE  LASER 


1.  INTRODUCTION 

The  C2n  Master  Database  was  assembled  in  support  of  the  Airborne  Laser  System 
Program  Office  (ABL  SPO).  The  database  consists  of  two  major  tables  of  data,  several 
forms  to  facilitate  data  entry  and  viewing,  some  queries  of  the  database,  and  some  reports 
based  on  the  database.  The  database  is  in  Microsoft  Access  format.  This  report  describes 
the  two  major  tables  in  the  database:  The  Master  Flights  Table  and  the  Altitude  Bins 
Table.  The  Master  Flights  Table  contains  the  date,  time,  and  location  of  each  balloon 
launch  and  various  parameters,  mostly  weather-related,  recorded  for  each  flight,  either 
determined  in  the  field  or  computed  later.  The  Altitude  Bins  Table  is  a  consolidation  of 
the  Thermosonde  optical  turbulence  data  from  each  flight.  A  third  table  entitled 
Campaign  is  a  handy  summary  of  each  campaign  and  the  status  of  the  data  reduction  for 
that  campaign.  This  report  describes  these  three  tables  in  detail,  including  the 
computational  methods  used  in  preparing  the  table  entries. 

2.  MASTER  FLIGHTS  TABLE 

By  selecting  the  Master  Flights  Table,  and  clicking  the  Open  button,  the  entire  table  can 
be  viewed.  To  view  the  data  for  a  single  flight,  open  the  Master  Flight  Viewer  Form. 

2.1.  Flight,  Date,  and  Time 

«A”  “B”  “C” 

Flight  Date  Time 

BAHSP001  6/1/97  8:18:00  PM 
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An  eight  digit  identifier  was  assigned  to  each  balloon  launch.  The  first  three 
characters  identify  the  location  where  the  launch  occurred;  the  next  two  characters 
indicate  the  season  during  which  the  launch  was  made;  the  sixth  digit  is  used  to  designate 
more  than  one  campaign  for  a  particular  season,  and  the  last  two  digits  number  each 

launch  made  during  the  particular  campaign. 

For  example:  BAHSP001  corresponds  to  Bahrain,  Spring,  launch  number  one; 

ROKWN102  is  Korea,  Winter,  launch  number  two;  RUHFA103  is  Saudi  Arabia,  Fall, 
launch  number  three;  and  HMNSP1 10  is  Holloman,  Spring,  launch  number  ten. 
ROKWN102,  RUHFA103  and  HMNSP1 10  all  designate  the  first  campaign  for  that 
Location.  This  “flight”  data  was  recorded  in  column  number  A.  The  “date”  of  each 


launch  was  recorded  in  column  B  and  the 

“time”  of  each  launch  (Local  time)  was  listed 

in  column  C. 

2.2.  ZF,  Winds,  Bottom  C2 ,  Top  C2 

“D”  “E” 

up> 

“G” 

ZF  Wind 

Bottom  C  2 

Top  C2 

30  y 

0 

30 

Column  D  contains  “ZF”  which  is  the  maximum  height  (in  kilometers)  that  the  balloon 


attained  in  the  launch  during  which  data  were  collected  for  C  \ .  A  value  of  0  km 

indicates  that  the  balloon  was  not  launched  successfully. 

If  wind  data  were  taken  it  is  shown  as  a  “y”  value  in  column  E;  if  no  wind  data  were 

taken  for  the  launch  it  is  shown  as  an  “n”  value  in  column  E. 

Bottom  C  l ,  in  column  F,  is  the  beginning  height  (km)  for  which  C  2  values  were 

taken  for  the  launch.  The  greatest  height  (km)  for  which  C  2  values  were  taken  is 

recorded  in  column  G. 
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2.3.  Trop(l),  Trop(2),  Trop(3) 

«jj»  cT  “Jp* 

Trop(l)  Trop(2)  Trop(3) 

17.5  0  0 

Columns  H,  I,  and  J  give  the  height  (km)  of  the  tropopause;  either  single,  double,  or 
triple.  The  tropopause  was  determined  by  examining  the  temperature  profile  of  the 
sounding.  Figure  1  shows  a  double  tropopause  for  the  Spring  98  campaign  at  Holloman 
AB,  New  Mexico.  As  depicted  by  the  temperature  curve  (solid  line)  there  is  an  abrupt 
change  in  the  decrease  in  temperature  at  13  km  and  again  at  19  km.  This  dramatic 
change  in  the  rate  of  temperature  decrease  gives  merit  to  a  Trop(l)  at  13  km  and  a 
Trop(2)  at  19  km. 


2.4.  Max  (m/s) 


“K” 


Max(m/s) 

33 


The  maximum  wind  speed  (m/s)  recorded  by  the  wind  data  taken  for  the  flight  is  given 
in  column  K. 

2.5.  200V(m/s),  200V  g  (m/s),  200V  v  (m/s),  200V  /V(m/s),  V  ag  Strength 

«L”  “M”  “N”  “O”  “P” 

200V(m/s)  200V  g  (m/s)  200V  ^  (m/s)  200Va^/V(m/s)  Vfl,  strength 
35  29.2  -5.8  -0.17  W 

The  next  five  columns  (L-P)  are  used  to  record  the  information  compiled  in 

determining  the  Ageostrophic  Intensity.  The  ageostrophic  wind  speed  was  computed  for 

all  thermosonde  launches  using  the  following  equation: 


V-V. 
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20 


80 


Relative  Humidity  (%) 

40  60 


*  |  1 1  1 1 1 1 1 1  *  | 1  •  j  1 1 1 1  j  1 1 1 1 1 1 1 1 1  j  1 1 1 1  j  |  j  | 

*50  -30  -10  10 

Temperature 


Temperature  (C)  vs.  Z  (Height  (Km))  Solid  Line  and  Relative  Humidity  (%) 
Z  (Height  (Km))  Dashed  Line  for  HMNSP107.  ' 


if  K  >  0,  Vg  >  V  and  where:  /  =  2Qsin  O  ;  Q  =  7.292<?^V  and  O  =  latitude. 

Where  V  is  the  observed  wind  speed,  Vg  is  the  geostrophic  wind  speed,  K  is  the  inverse 

of  the  radius  of  curvature  of  the  flow,  and  /  is  the  Coriolis  parameter  (20  sin  O)  K  is 
positive  if  the  flow  is  cyclonic,  negative  if  anticyclonic.  Equation  (1)  was  evaluated  for 
the  200  mb  pressure  surface  using  the  appropriate  latitude  of  each  thermosonde  launch 
when  evaluating  the  Coriolis  parameter.  V  is  taken  directly  from  the  NCEP  200  mb 
wind  speed  analysis.  The  radius  of  curvature  is  determined  graphically  from  the  NCEP 
200  mb  contour  analysis.  This  equation  was  evaluated  and  discussed  in  Newton  and 
Palmen  (1963)1  and  Palmen  and  Newton  (1969)2.  It  seems  suitable  for  our  initial 
approach  in  that  it  avoids  errors  associated  with  assuming  that  the  large-scale  analysis 
database  represents  an  accurate  measure  of  the  “observed”  wind  velocity  and  height 
contour  field  at  our  scales  of  interest.  The  curvature  factor  K  as  we  use  it  is 
applied  to  the  streamline  or  instantaneous  flow  field.  Strictly  speaking,  K  applies  to  the 
flow  trajectory.  The  difference  between  these  two  applications  at  200  mb  is  probably 
small  given  the  relatively  slow  movement  of  the  synoptic-scale  waves  relative  to  the 
observed  wind  speeds. 

The  ratio  of  (V  -  Vag)/V  was  then  calculated  for  each  profile  to  estimate  the  degree  of 

ageostrophic  forcing.  This  approach  is  patterned  after  the  diagnosis  of  Koch  and 
O’Handley  (1997)3  and  is  designed  to  be  a  crude  estimator  of  the  degree  of  unbalanced 
flow  associated  with  jet  stream  flow.  Qualitative  ageostrophic  forcing  categories  were 
assigned  to  this  ratio  as  follows:  the  range  0-0.2  is  assigned  as  weak,  0.2-0.4  is  moderate. 
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and  greater  than  0.4  represents  strong  ageostrophic  forcing.  Calculations  of  this  ratio 
show  the  highest  values  associated  with  the  strong  subtropical  jet  stream  or  jet  streak 
patterns.  The  magnitudes  and  sign  of  the  computed  ageostrophic  wind  speed  also  are 

consistent  with  those  found  in  Newton  and  Palmen  (1963). 

Figure  2.  is  the  200  mb  geopotential  height  chart,  from  the  NCEP  data  base,  used  to 
determine  the  radius  of  curvature  for  ROKSU022.  The  radius  (  R )  is  4.5  degrees  or, 
approximately,  495,000  meters  and  is  positive  due  to  the  cyclonic  curvature.  From  the 
NCEP  data  base,  the  200  mb  Speed  (  V  )  for  ROKSU022  is  10  m/s  and  the  coriolis 
parameter  (  /  )  for  Osan  AB,  Korea  is  0.000088  s'1  Substituting  these  values  in  equation 

(  1  )  after  solving  for  Vg  yields: 


V2 

v‘=7k+v 


Substituting  yields: 


(10  mis)2  .  1 

(,000088s-1  J  495000/W 


)  +  10(/w/s) 


Vg  =12.30  m/s 


(2) 


The  ageostrophic  speed  is  found  by  subtracting  the  geostrophic  speed  from  the  200 
mb.  speed  or: 

vn=v,-v  <3> 

Substituting  yields: 

Vag  =  12.30m/s-10.00w/s 

Vag  =  2.30m  Is 

The  intensity  of  the  ageostrophic  forcing  is  determined  by: 
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Intensity 


Substituting  yields: 


Intensity  = 


Intensity  =  .23 

Figure  3  is  the  200  mb  geopotential  height  chart  used  to  determine  the  radius  of 
curvature  for  ROKSU027.  In  this  case  the  radius  of  curvature  ( R  )  is  -4.5  degrees, 
negative  due  to  the  anticyclonic  flow.  From  the  NCEP  data  base,  the  200  mb  V  is  7  m/s 
for  ROKSU027  and  substituting  these  values  into  equation  (  2 )  yields  a  value  of: 

V 2 

v‘=Tk+v 

where:  V  -  V  =  -  (cyclonic);  and  V  —  V  =  +  (anticyclonic) 


V  =■  - L^  +  'jm/s 

g  (.000085-1)  495000 


Vg  =-1.12  +  7  =  5.88 


Subsequently,  substitution  into  equation  (  3  )  and  (  4  )  yields: 


V  =V  -V 

ag  '  g 


Vag  =  5.88-7.00 


Kg=- 1.12 


Intensity 


Intensity 


V 

-1.12 


Intensity  =  0.16 
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2.6.  Shear  Peaks 

“Q” 

Shear  Peaks 
10 

Column  Q  lists  the  number  of  vertical  shear  peaks  for  each  flight.  The  shear  peaks 
were  obtained  by  means  of  a  “wind  shear  calculator”  program.  The  program  uses  the 
“ascent  name.win”  files  from  each  balloon  ascent.  These  text  files  contain  three 
columns  of  data:  Altitude  (  km  ),  Wind  Speed  (  m/s ),  and  Wind  Direction  ( degrees  ). 
The  “wind  shear  calculator”  computes  and  constructs  graphs  of  the  wind  shear  with 
various  amounts  of  smoothing.  A  smoothing  of  the  wind  shear  data  points  is  necessary 
due  to  the  large  standard  deviation  between  individual  data  points.  Figure  4  is  a  plot  of 
vertical  wind  shear  results  from  four  ranges  of  smoothing  for  ROKFA1 15  (  Osan  AB, 
Korea,  14  Mar  1997 ).  It  was  decided  to  employ  an  1 1-point  smoothing  technique  for 
the  average  wind  shear  data  listed  here.  Figure  5  is  a  plot  of  the  1 1 -point  smoothing  of 
vertical  wind  shear  data  for  ROKFA1 15.  It  was  further  decided  to  take  values  of 
average  wind  shear  greater  than  .015  as  being  significant.  Column  Q  lists  the  number  of 
data  points  of  vertical  wind  shear  greater  than  .015  using  an  11-point  smoothing 
technique. 

2.7.  Jet  Speed  (  m/s  ),  Jet  Altitude  ( km  ) 

Jet  V  (m/s)  Jet  Alt  (km) 

38.3  14.1 

The  Jet  Speed  and  Jet  Altitude  were  recorded  from  the  wind  data  taken  during  each 
flight  and  are  recorded  in  columns  R  and  S.  The  value  of  Jet  V  differs  from  the  value  of 
max  (  m/s )  in  column  K  because  the  wind  speed  in  column  K  is  a  “smoothed”  value  of 
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0.04  - 


Figure  4.  Vertical  Shear  (/s)  vs  Altitude  (km)  at  7,  9, 11,  and  13  Point  Smoothing  for  ROKFA115. 
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Figure  5.  Vertical  Shear  (/s)  vs  Altitude  (km)  at  11  Point  Smoothing  for  ROKFA115. 


the  wind  speed  data  taken  during  the  flight.  Jet  V  and  jet  Alt  are  single  data  points  taken 
from  the  wind  profile  data.  Raw  data  points  were  used  to  get  the  best  determination  of 
where  the  balloon  was  located  as  it  entered  the  Jet  Stream  area.. 

2.8.  North,  South,  East,  West,  In  Jet,  Jet  Type,  Jet  (  m/s) 

«T»  U\V”  “Y” 

N  S  E  W  In  Jet  Jet  Type  Jet  (m/s) 

0  30  0  100  No  STJ  36 


Estimates  were  made  to  determine  the  position  of  the  balloon  launch  relative  to  the  Jet 
Stream.  A  determination  of  where  the  balloon  was,  in  respect  to  where  it  was  launched, 
as  it  reached  the  Jet  Stream  altitude  ( column  S  )  was  made  using  the  “wind  shear 
calculator”  program  discussed  earlier.  Figure  6  is  a  plot  of  the  distance  traveled  for 
ROKFA1 15  as  generated  by  the  “wind  shear  calculator”  program  and  shows  the 
trajectory  of  the  balloon  after  launch.  The  NCEP  data  base  was  used  to  obtain  plots  of 
200  mb  winds  for  each  flight  (  see  Figure  7  ).  The  200  mb  wind  plots  were  used  to 
determine  the  location  of  the  jet  stream  and  the  computed  position  of  the  balloon  was 
added  to  the  plot.  The  distance  from  the  balloon  to  the  core  of  the  Jet  Stream  was 
measured  (  assuming  110  km  distance  for  each  degree  lat/long  )  and  recorded  in 
columns  T  -  W  as  the  distance,  in  kilometers,  N,  S,  E,  W  of  the  Jet  core.  Column  X  is 
used  to  record  whether  the  balloon  is  in  the  Jet  Stream  or  not.  For  a  Polar  Jet,  a  200  mb 
wind  of  >  40  m/s  was  used  to  determine  if  the  balloon  was  in  the  Jet  Stream,  while,  for  a 
Sub-Tropical  Jet,  a  200  mb  wind  of  >  25  m/s  was  the  determining  factor.  The  200  mb 
wind  speed  is  recorded  in  column  Z  (  jet  (  m/s  )).  The  200  mb  wind  speed  was  taken 
from  the  NCEP  data  base  (  see  Figure  7  ).  Finally,  the  Jet  Type  was  determined  by  the 
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.  200  mb  Winds  (m/s)  for  R0KFA1 15  Showing  Site  Launch  and  Position 
of  Balloon  at  200  mb.  Level. 


location  of  the  Jet  Stream,  as  shown  on  the  NCEP  200  mb  wind  plot,  and  recorded  in 
column  Y.  If  the  Jet  Stream  was  located  North  of  40  degrees  latitude,  it  was  classified  as 
a  Polar  Jet  and  if  it  was  located  South  of  40  degrees  latitude  (  Northern  Hemisphere  ),  it 
was  classified  as  a  Sub-Tropical  Jet.  There  were  cases  when  the  two  Jets  combined  and 
this  was  also  recorded. 

2.9.  Weather,  Surface  Temperature  (  C ),  Frontal  Passage,  Sea-Level  Pressure  (  mb  ), 
Wind  (  m/s  ),  Relative  Humidity  (  %  ) 

“AA”  “AB”  “AC”  “AD”  “AE”  “AF” 
Weather  SfcT(C)  FROPAS  SLP  ( mb )  wind  ( m/s )  RH ( % ) 
Partly  Cloudy  39  No  1004  2  15 

Columns  AA  to  AF  contain  the  synoptic  weather  conditions  for  each  flight  as  close  to 

the  time  of  launch  as  possible.  The  data  was  obtained  mostly  from  the  NCEP  data  base. 

Surface  temperature,  sea-level  pressure,  surface  wind  speed  and  relative  humidity  were 

taken  directly  from  the  NCEP  data  base  while  the  weather  and  frontal  passage  data 

were  either  inferred  from  the  NCEP  data  base  or  acquired  from  the  synoptic  notes  taken 

for  about  half  of  the  flights  by  a  meteorologist  at  the  site.  These  notes  are  also  recorded 

in  the  comments  contained  in  the  last  column  (  AO ). 

2.10.  Jet  Enhancement,  Tropospheric  Enhancement,  Stratospheric  Enhancement 

“AG”  “AH’  “AT’ 

Jet  enhance  trop  enhance  strat  enhance 
no  strong  no 

A  determination  was  made  as  to  whether  there  was  enhancement  of  C  \  caused  by 

the  Jet  Stream,  the  Tropopause  or  the  Stratosphere.  If  there  was  enhancement,  it  was 
recorded  as  strong  or  weak  in  columns  AG,  AH,  and  AI.  The  enhancement  was 
determined  by  comparing  the  plot  of  C^  vs  Z  to  the  CLEAR  1  model.  Figure  8  is  the 

plot  of  C  l  vs  Z  for  flight  HMNSP104  ( 14  May  1998  )  taken  at  Holloman  AB,  New 
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Mexico.  For  this  flight,  there  is  no  evidence  of  jet  enhancement  ( taken  from  column  X 
which  has  a  negative  since  200  mb  wind  speed  was  ^  25m/s),  but  the  values  of  C  2  are 
significantly  larger  at  the  tropopause  region  (13  and  19  km.  from  columns  H  and  I )  than 
the  CLEAR  1  model.  Above  the  tropopause  (  20  km.  ),  values  for  C  2  are  lower  than 

CLEAR  1,  therefore,  there  is  no  stratospheric  enhancement. 

2.11.  Deformation,  Intensity 

“AT  “AK” 

deformation  intensity 

4.37  Light 

A  deformation  term  was  calculated  for  each  flight  and  is  shown  in  column  AJ.  The 
deformatiom  term  was  used  to  obtain  a  measure  of  the  horizontal  shear.  From  Saucier 
(  1955  )4  the  following  equation  was  used: 


\DEF\  = 


rdu_ 

Kdx 


dv 

fy. 


dv  du ^ 
{dx  +  dy^ 


(5) 


Appropriate  values  of  u  and  v  were  obtained  from  the  NCEP  data  base.  Figure  9  is  the 
plot  of  u  for  ROKSP101  (  Osan  AB,  Korea,  5  May  1998 ).  A  circle  of  radius  2  degrees 
was  drawn  at  the  position  of  the  launch  and  values  of  the  change  of  u  in  the  x  (  W-E  )  and 
y  (  S-N  )  direction  was  measured  from  the  plot.  For  flight  ROKSP 101,  du/dx  is  40  -  28 
or  12  m/s  and  du/dy  is  28  -  40  or  -12  m/s.  In  a  similar  manner  values  of  v  were  obtained 
(  see  Figure  10 ).  For  ROKSP101,  dv/dx  is  12.5  -  1  or  1 1.5  m/s  and  dv/dy  is  -1.25  - 
14.5  or  -15.75.  Substituting  these  values  into  equation  ( 5  )  yields: 


\DEF\  =  [(l  2  -  (- 1 5.75))2  +  (l  1 .5  + 

| Off]  =  [(27.75)!  +  (-  5)’]!  and,  \DEF\  =  27.75. 
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Based  on  results  from  Mancuso  and  Endlich  (  1966  )5 ,  values  up  to  10  were  considered 
light  deformation,  between  10  and  20  were  moderate  and  greater  than  20  were  considered 
strong  deformation.  The  intensity  of  deformation  is  listed  in  column  AK. 

2.12.  Mountain  Wave  Intensity 

“AL” 

Mtn  wave  inten 
Light 

Column  AL  lists  the  mountain  wave  intensity  for  flights  where  mountain  ranges  were 
present  and  the  possibility  of  a  mountain  wave  existed.  For  this  data  base,  mountain 
waves  were  examined  for  Korea  and  Holloman.  The  approach  employed  was  to  use  a 
nomogram  from  the  Air  Force  Weather  Agency  manual  AFWA  TN-98/0002,  15  Jul 
1998,  Meteorological  Techniques  (  see  Figure  11 ). 

The  mountain  wave  nomogram  method  predicts  mountain  wave  intensity  qualitatively 
in  categories:  none,  light,  moderate,  and  severe  as  a  function  of  maximum  wind  speed 
below  500  millibars  and  sea  level  pressure  differential  across  the  mountains.  For  our 
analysis,  the  pressure  differential  is  calculated  simply  by  subtracting  the  east  Korean 
coast  sea  level  pressure  from  the  west  Korean  coast  sea  level  pressure  determined  from 
weather  maps.  The  500  mb  and  below  wind  data  are  taken  from  thermosonde 
measurements.  For  Holloman,  the  boundary  (  E-W )  of  the  state  of  New  Mexico  was 
used  to  determine  pressure  differential.  The  wind  direction  was  also  taken  into  account 
since  a  wind  perpendicular  to  the  mountain  range  would  give  a  full  effect,  while  a  wind 
parallel  to  the  mountain  range  would  have  no  effect.  For  Korea  and  Holloman  a 
westerly  wind  (  270  degrees )  was  considered  to  be  perpendicular  to  the  mountain  range. 
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MOUNTAIN-WAVE  NOMOGRAM 


PRESSURE  (mb) 

KNOTS  s  Maximum  Wind  Speed  Beiow  500mb 

PRESSURE  =  Sea-level  Pressure  Differential  Across  Mountains 


Figure  11.  Mountain  -  Wave  Nomogram  Used  to  Determine  Mountain  -  Wave  Intensity. 
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2.13.  Aircraft,  RFI,  Comments 


“AM” 

“AN” 

“AO” 

aircraft 

RFI 

comments 

n 

n 

(  as  stated) 

If  an  aircraft  was  employed  to  supplement  the  flight  it  is  entered  in  column  AM.  In 
some  launches,  radio  frequency  interference  (  RFI )  was  a  problem  and  if  this  was  the 
case  it  was  noted  in  column  AN.  The  last  column  (  AO  )  contains  comments  made  by 


the  launch  team  at  the  time  of  the  launch.  These  comments  have  been  preserved  from 
their  initial  entry. 

3.  OPTICAL  TURBULENCE  DATA  REDUCTION  AND  STORAGE 
Validated  data  from  individual  balloon  ascents  has  been  averaged  into  altitude  bins  and 
stored  in  the  table  “Altitude  Bins”  in  the  Master  Database. 

3.1  Data  Source 

The  source  of  the  turbulence  data  for  the  Database  is  the  standard  “dot-text” 
files  for  the  individual  ascent  profiles.  This  file  is  not  the  raw  data  from  the  radiosonde, 
but  one  of  the  many  files  created  in  the  data  reduction  process.  This  particular  file  is  in  a 
form  that  is  usually  required  for  atmospheric  characterization  for  optical  propagation 
calculations.  The  files  are  validated  in  the  sense  that  they  are  carefully  checked  and 
known  file  defects  such  as  excessive  noise  and  baseline  drift  are  eliminated.  An 
excerpt  from  a  text  file  follows: 


0.0116 

1022.8998 

3.8000 

68.0000 

4.8134e-015 

0.0227 

1021.5002 

3.8000 

68.0000 

1.2535e-014 

0.0323 

1020.2997 

4.1000 

68.0000 

1.6684e-014 

0.0379 

1019.5999 

4.3000 

68.0000 

1.0248e-014 

0.0451 

1018.7002 

4.6000 

67.0000 

5.6553e-015 

0.0515 

1017.9002 

4.8000 

67.0000 

3.7696e-015 
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The  columns  are:  altitude  (km)  above  sea  level ,  pressure  (mb),  temperature  (°C), 
relative  humidity  (%),  and  the  structure  constant  for  index  of  refraction,  C„2  (m 2  ).  The 
lines  represent  data  for  each  line  of  telemetry,  which  is  received  at  intervals  of  1  to  2 
seconds,  depending  on  the  radiosonde  used  for  the  experiment.  These  files  can  vary  in 
length  from  3000  to  6000  lines,  and  contain  much  more  precision  than  is  usually  required 
for  optical  calculations.  For  any  ascent,  there  may  be  missing  values  at  the  beginning,  at 
the  end,  or  anywhere  in-between.  This  could  be  due  to  telemetry  dropouts  or  other  problems. 
3.2  Data  Peparation  for  Table 

One  of  the  first  decisions  that  had  to  be  made  to  assemble  the  database  was  the  best 
thickness  in  altitude  to  be  used  for  storage.  By  averaging  the  Q,2  for  larger  and  larger 
altitude  bins,  it  was  determined  that  the  precision  of  a  Rytov  variance  calculation  did  not 
begin  to  seriously  degrade  until  the  bins  became  thicker  than  500m.  After  further 
consideration,  it  was  decided  to  use  a  150m  bin  size,  which  is  the  same  size  that  is  used 
for  WSMR  radar  data.  Further,  the  bins  were  all  aligned  such  that  identical  altitude 
ranges  would  be  used  for  each  bin,  so  that  such  calculations  such  as  the  mean  turbulence 
for  a  given  altitude  could  be  easily  computed.  Finally,  the  altitudes  used  were  chosen  to 
be  identical  to  those  used  for  the  WSMR  radar  data.  The  bottom  altitude  for  the  first  bin 
is  3.222  km,  the  top  of  the  first  bin  is  3.372  km,  and  the  nominal  altitude  in  the  file  is  the 
mean,  3.297  km.  The  last  bin  was  at  a  nominal  altitude  of  29.997  km,  since  balloon  data 
were  typically  disregarded  above  30  km.  The  balloons  usually  burst  around  that  altitude, 
and  there  is  a  question  about  the  accuracy  of  the  measurements  above  that  altitude  due  to 
heat  transfer  considerations. 

The  Thermosonde  actually  measures  the  structure  constant  for  temperature,  CT2.  This 
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quantity  is  then  converted  to  C„2  using  the  formula: 

Cl=[l9xW*j^Cl  (6) 

where  P  is  pressure  in  mb,  T  is  temperature  in  K,  and  Cj2  is  in  K2m"2/3 

The  atmospheric  conditions  in  the  multiplier  impose  decreasing  Cn2  values  with 

altitude,  while  minimum  and  maximum  Or2  excursions  remain  relatively  constant  with 

altitude.  Therefore  it  was  decided  to  convert  Q,2  back  to  Or2  prior  to  averaging  over 
altitude.  The  average  for  bin  j  is  the  mean  value  of  an  assumed  continuous  function 

defined  by  the  formula: 


J  C}(z,)dz 


2  _  zfl 


c2  = 


Z-  —Z; 
J 1 


where  z  is  the  altitude,  Zji  is  the  minimum  altitude  within  the  bin  j,  and  Zj„  is  the 
maximum  altitude  in  the  bin  j.  The  numerical  integration  is  performed  using  the 
trapezoid  rule.  Since  the  lower  and  upper  altitudes  of  the  bin  do  not  usually  occur  at  the 
data  altitudes,  the  values  for  the  Zji  and  Zj„  are  estimated  by  linear  extrapolation  of 
the  bounding  lines  of  data.  There  must  be  at  least  three  lines  of  data  within  the  bin  for  the 
bin  to  be  considered  in  the  database.  If  a  bin  is  not  valid  for  lack  of  the  required  data, 


the  bin  is  not  entered  into  the  database. 

Means  are  also  calculated  for  temperature  and  pressure,  with  the  purpose  of 
recomputing  a  valid  value  of  C„2.  With  this  in  mind,  it  is  the  mean  value  based  on  P2  and 
the  mean  value  based  on  T4  that  is  computed.  Each  line  in  the  Altitude  Bins  table, 
contains  the  balloon  flight  identifier,  and  mean  values  of  altitude,  Ct  ,  P,  T,  and  Cn  . 


Note  that  the  temperature  mean  (T-4)  is  computed  using  degrees  Kelvin,  then  it  is 
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converted  back  to  Celsius.  There  is  also  a  unique  %  number  automatically  assigned  to 
each  entry  line  by  the  program.  A  sample  of  the  table  is  shown  below. 


Table  2.  Table  of  Flight  Bins  and  Computed  Means. 


H  1 1 

liWUlii 

iSlI 

IfifggMTOTfiil 

mwzmm 

ilHUH 

64433 

roksu020 

12.897 

0.0000163820 

187.1999; 

-50.2838: 

1 .4523E-1 8 

64434 

roksu020  ! 

13.047 

0.0000141660 

182.9569: 

-51.4488: 

1 .22491-18 

64435 

roksu020  | 

13.197 

0.0000098826 

178.78951 

-52.4193; 

8.30481-19 

64436 

roksu020 

13.347 

6.66661 i 0490 

174.6991; 

-53.6162; 

9.0608E-19 

64437 

roksu020  1 

13.497 

6.6666689451 

170.67891 

-54.7668; 

7.1503E-19 

:  64438 

roksu020 

13.647 

6.6666673323 

166.7330; 

-55.8024! 

5.7006E-19 

In  a  typical  use  of  the  database,  a  selection  is  made  for  either  a  single  flight,  or  a  group 
of  flights  based  on  conditions  in  the  Master  Flights  Table.  Next  the  desired  data  are 
extracted  from  the  Altitude  Bins  Table  using  a  query.  Finally  calculations  are  made 
using  the  extracted  data.  A  simple  example  calculation  is  the  mean  Cn  for  each  altitude 
bin,  which  can  be  accomplished  using  the  AVG  function  built  into  the  database. 

4.  CAMPAIGN  SUMMARY  TABLE 

The  campaign  summary  table  consists  of  the  following:  campaign  designator  for  the 
flight  identification  number,  campaign  location,  date  of  first  and  last  balloon  launch  of 
the  campaign,  and  number  of  flights  entered  into  the  database.  There  are  also  check 
boxes  to  indicate  completion  of  data  entry  into  the  Master  Flight  Table,  the  Altitude 
Bins  Table,  if  the  campaign  had  contemporaneous  airborne  data  collections,  and  if  the  data 
had  been  through  the  quality  control  procedure.  A  convenient  way  to  view  these  data  is 
to  call  up  the  Campaign  Summary  Form,  which  displays  the  data  from  each  campaign  on 
a  single  window. 
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5.  CONCLUSION 


In  compiling  the  Master  Data  Base,  every  attempt  was  made  to  use  thermosonde 
measured  values  wherever  possible.  In  doing  this,  it  is  hoped  that  the  accuracy  of  the 
Master  Data  Base  could  be  kept  at  the  highest  possible  level.  However,  some  of  the 
parameters,  i.e.,  deformation,  ageostrophic  wind,  mountain  wave, etc.,  had  to  be 
computed  for  large  scale  analyses.  Where  this  was  necessary,  methodologies  were 
employed  that  retained  the  highest  possible  accuracy  and  the  results  were  considered 
as  accurately  as  possible.  For  example,  actual  deformation  values  were  determined  but 
were  then  recorded  as  qualitative  values  as  light,  moderate  and  strong  intensities. 
Furthermore,  all  the  parameters  that  could  be  garnered  from  the  actual  flight 
data  were  used  and  recorded  in  raw  form.  This  allows  future  analysis  to  be 
performed  should  new  ideas  evolve  in  the  overall  ABL  study.  The  Master  Data  Base  is 
considered  a  good  summary  of  the  numerous  weather  related  parameters  occurring  during 
the  entire  ABL  field  program  to  date  and  can  easily  be  updated  to  include  additional 
research  data.  The  Altitude  Bins  Table  contains  optical  turbulence  parameters  Q2  and 
Ct2  along  with  the  mean  altitude,  mean  pressure  and  mean  temperature  of  each  flight  bin 
with  at  least  3  interior  data  records.  The  Campaign  Summary  Table  is  a  good  first  table 
to  use  to  determine  what  campaigns  are  part  of  the  database,  and  the  status  of  the  data 
reduction.  The  C„2  Master  database  also  contains  many  forms  that  are  reports  which  should 
be  useful  in  viewing  and  reporting  information. 
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